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FINITE  ELEMENT  METHODS  FOR  THE  SOLUTION 
OF  PROBLEMS  WITH  ROUGH  INPUT  DATA 

V  -f  # 

I.  Babuska  and  J.E.  Osborn 

1.  Introduction 

The  finite  element  method  is  based  on  the  variational  or  weak 
formulation  of  the  boundary  value  problem  under  consideration.  The 
aporoximate  solution  is  obtained  by  restriction  of  the  variational 
formulation  to  finite  dimensional  trial  and  test  spaces. 

Accuracy  of  the  approximate  solution  is  achieved  by  the  choice  of 
a  trial  space  with  good  approximation  properties  and  the  choice  of  a 
test  space  which  guarantees  that  the  finite  element  solution  is  an 
approximation  of  the  same  quality  as  is  the  best  approximation  of  the 
exact  solution  by  the  trial  functions. 

The  success  of  the  finite  element  method  as  a  practical  coiqputational 
tool  is  related  to  the  special  construction  of  the  trial  and  test  func¬ 
tions  in  terms  of  "element"  trial  and  test  functions  (defined  on  the 
finite  elements)  satisfying  appropriate  constraints  at  the  nodes.  This 
element  based  structure  is  the  basis  for  the  architecture  of  existing 
finite  element  codes.  Usually  the  "element"  trial  and  test  functions 
are  polynomials. 

It  is  well  known  that  when  the  solution  is  rough,  which  will  in 
aeneral  be  the  case  when  the  data  is  rough,  then  the  use  of  polynomials 
as  "element"  trial  functions  does  not  lead  to  good  accuracy.  In  this 
paper  we  will  construct  "element"  trial  functions  which  reflect  the 
properties  of  the  problem  under  consideration.  This  will  lead  to 
Increased  accuracy  -  the  maximum  possible  accuracy  -  while  not  changing 
the  structure  of  the  code.  In  fact,  the  approach  can  be  implemented 
by  changinn  only  the  "element"  trial  functions  employed  in  the  method. 

The  questions  of  optimal  trial  functions  relative  to  a  class  of  problems 
is  also  addressed. 

The  abstract  framework  of  the  approach  is  given  in  Section  2.  In 
Section  3,  five  examples  of  various  types  are  presented  and  in  terms 
of  these  examples  the  ideas  presented  in  Section  2  are  developed.  In 
Section  4,  we  make  some  general  comments  on  the  design  of  finite 
element  methods  for  problems  with  rough  coefficients. 

Detailed  proofs  of  the  results  in  the  paper  will  be  presented 
elsewhere.  The  ideas  in  the  paper  can  be  generalized  in  various 
directions. 
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Throughout  the  paper  we  will  use  the  usual  spaces  and  Sobolev 

spaces 


2.  Variational  Approximation  Methods 

In  this  section  we  discuss  variational  approximation  methods  and 
state  a  basic  error  estimate. 

Let  and  be  two  Hilbert  spaces  and  let 

'^l,h''l*''Hj^  0  <  h  <  1,  H2,h' 

families  of  Hilbert  spaces  satisfying 


'H 


2,h 


,  0  <  h  <  1,  be  two  one  parameter 


^1  ”l,h' 


l,h 


and 


"2  =  ''2,h'"“"H 


2,h' 


llullu  Vh  and  Vu  €  H, 
«!  1 


lulU  Vh  and  Vu  €  H- . 

”2  2 


(2.1) 


(2.2) 


Let  Bjj(u,v)  be  a  bilinear  form  on  ^1  h  ^2  h  satisfying 

|b.  (u,v)|  <  M))u|)^  )|v)|^  Vu  €  h  /  and  Vh 

l,h  2,h  '  ' 


(2.3) 


(M  is  independent  of  h) . 

Let  f  be  a  bounded  linear  functional  on  Corresponding  to 

f  we  assume  that  there  is  a  unique  Uq  €  satisfying 


Bh(Uo»v)  =  f(v)  Vv  €  and  Vh. 


(2.4) 


Uq  can  be  thought  of  as  the  exact  solution  to  a  boundary  value  problem 
under  consideration  and  is  unknown.  Bj^  and  ^  are  given  input 

date  and  are  known. 

We  are  interested  in  approximating  u^  and  toward  this  end  we 
assume  we  are  given  or  have  chosen  finite  dimensional  spaces  c  ^ 

and  ^2  j,  with  dim  Sj^  =  dim  S2  ^  and 


inf 


“«®l,h  “'®2,h 


'lull 


:iv 


sup  I 

=1 


B.  (u,v)  I  >  Y  >  0  Vh 


(2.5) 


H 


2,h 


B-.- 


(y  is  independent  of  h)  .  Then  we  define  Uj^  by 


Di'trlbutlon/ 


.« vaiiablllty  Cod 
[Avail  and/oa 


Di;it 


Special 
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^4 


m 


y 

I  v't 


'Vi 

■,'.S 
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I 

■'?vf 


%  ^  ®l,h 


Bh(Uh,v)  =  Vv  ^ 


(2.6) 


and  consider  u.  to  be  an  approximation  to  u».  (2.6)  is  uniquely 


solvable.  We  often  denote  u^  be  u.  Note  that  u^^  is  defined  for 


every  ^  h‘  Given  bases  for  Sj^  and  S2  (2.6)  is  reduced 
to  a  system  of  linear  equations.  But  (2.6)  does  not  define  an  algorithm 
for  finding  Uj^  since  it  depends  on  the  unknown  Uq.  We  note  from 
(2.4),  however,  that  Bjj(Uq,v)  =  f(v)Vv  €  S2  provided 


S2,h  ^  «2 


(2.7) 


We  now  assume  that  (2.7)  holds.  (2.6)  can  then  be  written 


'^h  «  ^l,h 


Bh(Uh,v)  »  f(v)  Vv  €  S2 


(2.8) 


(2.8)  is  called  a  variational  approximation  method. 

Having  defined  Uj^  we  are  interested  in  an  estimate  for 

This  is  provided  by  the  standard 


Theorem.  The  error  satisfies 


'"0  -  "h"H. 


<  (1  +  Y  M)  inf  l|u-  -  xl 


(2.9) 


where  M  and  y  are  the  constants  in  (2.3)  and  (2.5),  respectively. 


For  a  complete  discussion  of  variational  approximation  methods  see  [1,2] 
The  spaces  Sj^  are  called  trial  spaces  and  the  functions  in 
are  called  trial  or  approximating  functions.  The  spaces  S2 


are  called  test  spaces  and  the  functions  in  S.  :  are  called  test 

2,h 


functions . 


Remarks:!)  In  our  applications  there  will  be  a  bilinear  form  B(u,v) 
defined  on  and  satisfying  Bj^(u,v)  =  B(u,v)  Vu  €  v  (  H^, 


Prom  (2.4)  we  see  that 


B(Uq,v)  =  f(v)  Vv  ( 


(2.10) 


(2.10)  is  the  variational  formulation  of  our  boundary  value  problem  and 
Hj^,  H2,  and  B  are  the  spaces  and  the  bilinear  form  in  this  formulation. 
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2)  (2.9)  suggests  we  choose  S,  .  so  that  inf  ||Uq  -  x''u  is 

' 

small,  i.e.,  so  that  the  trial  functions  have  good  approximation  prop¬ 
erties,  and,  with  so  chosen,  S2  can  be  selected  so  that  (2.5) 

holds  with  as  large  a  constant  y  as  possible. 

3)  In  many  applications  we  can  choose  h  ^1'  others  the 

requirement  that  inf  l!u^  -  xHu  small  leads  one  to  choose 

x€Si,h  ”l,h 

Si  h  4  The  trial  space  is  then  nonconforming  in  the  sense  that 

Si  h  basic  variational  space  Hj .  This  fact  leads 

to  a  use  of  the  family  of  spaces  ^  and  forms 

In  certain  situations  one  has  and  one  wishes  to  choose 

Sj  h  “  ®1  h’  i^  h  nonconforming,  S2  will  be  also, 

and  we  are  led  to  the  use  of  the  family  Note  that  in  this 

circumstance, 

Bh(Uo»v)  =  f(v)  Vv  €  S2  (2.11) 


is  not  valid.  Nonetheless,  an  approximation  Uj^  can  still  be  defined 
as  in  (2.8).  This,  in  fact,  is  what  is  done  in  the  class  of  methods 
commonly  referred  to  as  nonconforming  in  the  finite  element  literature 
(see,  e.g..  Section  4.2  in  [4]).  The  error  analysis  for  such  problems 
does  not  follow  directly  from  (2.9)  and  the  additional  conqplications 
in  the  analysis  are  due  to  the  fact  that  (2.11)  does  not  hold. 

In  the  methods  discussed  in  this  paper  we  will  always  have 
^2  h  ^  ^2'  test  spaces  will  be  conforming.  We  emphasize 

that  the  choice  of  nonconforming  trial  spaces  causes  no  difficulty  in 
the  analysis,  provided  the  test  spaces  are  conforming. 


3.  Examples 

In  the  section  we  will  discuss  the  approximate  solution  of  five 
specific  boundary  value  problems  with  rough  input  data.  In  each 
examole  we  will  use  a  variational  approximation  method  employing  trial 
functions  which  reflect  the  properties  of  the  underlying  problem  in 
that  they  provide  an  accurate  approximation  to  the  unknown  solution. 

a .  A  Two  Point  Boundary  Value  Problem  with  a  Rough  Coefficient 
Consider  the  problem 

Lu-  -  -(a(x)u' ) '  *  f,  0  '  X  "  1 

°  (3.1) 

Uq(0)  -  Uq(1)  =  0 
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where  a(x)  is  a  rather  arbitrary  function  satisfying  0<o<a(x)  ^  R 
and  f  €  L2(0,1).  This  simple  model  problem  arises  in  the  analysis 
of  the  displacements  in  a  tapered  elastic  bar.  f  represents  the  load, 
a(x)  the  elastic  and  geometric  properties  of  the  bar,  and  Uq  the 
displacement.  If  the  bar  has  smoothly  but  rapidly  or  abruptly  varying 
(as  in  the  case  of  composite  materials)  material  properties,  then 
a(x)  will  be  a  smoothly  but  rapidly  or  abruptly  varying  function,  i.e., 
a  rough  function. 

The  variational  formulation  of  (3.1)  is 


Ug  €  ftj^(0,l) 

*6 


f  auiv'dx  =  [  fvdx  Vv  €  H. (0,1). 

J  0  Jo 


(3.2) 


This  has  the  form  (2.10)  with  ~  ^1'  ®(“»v)  =  au'v'dx,  and 

fdx. 

J  0 


f  (V) 


It  is  Icnown  that  the  standard  finite  element  method 


enploying  c”  piecewise  linear  trial  and  test  functions  (i.e.,  the 
variational  approximation  method  determined  by  the  form  B  and 


piecewise  linear  trial  and  test  spaces)  does  not  yield  accurate  approxi¬ 
mations  to  the  solution  of  (3.2)  (or  (3.1))  when  a(x)  is  rough.  We 
thus  consider  an  alternate  method. 

Let 


{0  =  Xg  <  Xj^  <. .  .<  x^  =  1}  be  an  arbitrary  mesh  on 


h  “  0  ‘'I 

[0,1]  and  set  Ij  =  “  *j"*j-l'  ^ 


points  Xj  are  called  nodes.  For  the  trial  space  we  choose 


S,  u  =  {os  For  each  j,  u 
l,h 


-  a  linear  combination  of  1  and 


dt  . 
iTtT' 


is  a  solution  of  (au')'  ^  0,  (3.3) 

u  is  continuous  at  the  nodes,  and  u(0)  »  u(l)  »  0} 


I . 
3 


and  for  the  test  space 

S2  *  {v:  For  each  j  ,  v 

u  is  continuous  at  the  nodes,  and  u(0)  =  u(l)  =  0}.  (3.4) 


=  a  linear  combination  of  1  and  x. 


I . 
J 


Now 


h  '■  ^1  ®2  h  ^2  thus  choose 


l,h 


”1'  ^2,h  “  ^2'  ®h  “ 


(3.5) 


However,  we  could  also  choose 


=  ^f2,h  ~  ^l,h  ~  ^  H^(Ij)  for  each  j,  u(0)  =  u(l)  =  0}, 

l|ulL  =  !|u!l„  =  l|ul|«  =  [[  u^dx  +  I  [  (u')^dx]^^^,  (3.6) 

”l,h  ”2,h  “l,h  Jq  j=l 

Bj^(u,v)  =  Ij  I  u'v'dx. 


With  either  the  choice  (3.5)  or  (3.6),  we  see  that  (3.2)  has  the  form 
(2.4)  and  that  (2.3)  is  satisfied  with  M  =  B.  We  will  use  the 

choice  (3.6)  since  it  leads  more  naturally  to  the  choice  dictated  to 
us  in  Example  b  below. 

We  next  consider  the  variational  approximation  method  (2.8)  deter¬ 
mined  by  ®h'  ®1  h'  ®2  h' 


"h  ^  ®l,h 


(Uh,v)  .  r 


au^v’dx  =  j  fvdx  Vv  €  Sj  jj. 


(3.7) 


Regarding  (3.7)  one  can  prove  the  inf-sup  condition  (2.5),  namely 

inf  sup  I  I  au'v'dxl  =  Y(a,B)  >  0  Vh 

u€S,  V,  v€S-  V.  “'0 

2,h  (3.8) 


H,=i  iluHfl 


Hi=l 


and  the  approximation  result 


inf  llug  -  xHg  ”5  N(a,0)hl|(au‘)  MIl  =  N(a,B)h||f  l|j^^ 
[€S,  .1  22 


(3.9) 


where  y{ci,Q)  and  N(a,B)  depend  on  a  and  B  but  not  otherwise  on 
a(x)  nor  on  h.  From  (3.8),  (3.9)  and  the  basic  estimate  (2.9)  we 
immediately  have 


Theorem  a.  The  error  u,.  -  u,.  satisfies 
-  0  h 

lu^  -  <  C(a  ,e)hi|fi!j^^ 


(3.10) 


where  C(oi,6)  =  (1  +  y  ^  (a, 3) M(a,  3)  ]N (a,  6)  depends  only  on  a  and  6. 
Thus  we  have  a  1—  order  estimate  for  ^'Uq  -  which  is  uniform 


over  all  a(x)  satisfying  a  5  a(x)  <  3.  For  a  proof  of  (3.8)  and 

(3.9)  see  [3].  The  estimate  (3.10)  is  the  best  possible  estimate 
for  f  €  L2  This  can  be  seen  by  an  application  of  the  theory  of 
N-widths  (see,  e.g.,  [2,  8]). 

(3.9)  shows  that  as  defined  in  (3.3),  yields  accurate 

approximations  to  the  solutions  u^  of  (3.1).  (3.8)  shows  that  (2.5) 

holds  with  our  choice  of  and  S2  as  defined  in  (3.4).  Thus 

our  trial  and  test  spaces  have  been  chosen  in  accordance  with  the 
suggestions  in  Remarlc  2  in  Section  2.  We  further  note  that  if 
<t>i  , . . .  usual  basis  for  S~  .  defined  by  (x. )  =  6..  and 

IN  ^  f  n  ^ijij 

, . . .  is  the  basis  for  S,  ,  defined  by  (x. )  =  6..,  then  the 

iN  ijij 

matrix  of  (3.7)  (the  stiffness  matrix)  is  given  by 


I  a1i'^((ij^dx  =  I  aj^((i^<|»|dx 


(3.11) 


where  a^  is  the  piecewise  harmonic  average  of  a(x) ,  i.e.,  the  step 
function  defined  by 


Thus  the  stiffness  matrix  is  synmetric  (which  is  not  immediate  since 
we  are  using  different  trial  and  test  spaces)  and  is  as  easily  computed 
as  is  the  stiffness  matrix  for  the  standard  method  employing  S2  for 

both  trial  and  test  space,  namely  j  ajj())^(J> j^dx,  where  a^^  is  the  piece- 

0 

wise  average  of  a(x). 

The  accuracy  and  robustness  of  the  approximation  u^  is  further 
shown  by  the  estimate 

|u(x.)  -u.(x.)|  <  C(a,8)vJ(a)h^llfll  ,  j=l,...,n-l  (3.13) 

J  "  J 

where  VQ(a)  denotes  the  total  variation  of  a(x).  Thus  we  have 
second  order  convergence  at  the  nodes  even  if  a(x)  has  several  jumps. 
The  proof  of  (3.13),  which  does  not  follow  the  lines  suggested  by 

(2.9),  can  be  found  in  [3]. 


b.  A  Special  Class  of  Two  Dimensional  Boundary  Value  Problems  with 
Rough  Coefficients 
Consider  the  problem 

Lu  =  -(a(x)u  )  -  (b(y)u  )  =  f(x,y),(x,y)  €  Q  =  [0,1]^ 

XX  y  y 

u  =  0  on  r  =  (3.1^) 

where  0  <  a  <  a(x) ,b(y)  <  3  and  f  €  L2.  This  problem  generalizes 
the  model  problem  considered  in  a  .  The  variational  formulation  of 
(3.14)  is 

u  € 

B(u,v)  =  f(v)  \/v  €  Hj^(n)  (3.15) 


where 

B(u,v)  =  I  (au^v^  +  bUyVy)dxdy 
fi 

and 

f(v)  =  I  fvdxdy. 

n 


This  is  of  the  form  (2.10)  with 

Let  be  a  uniform  triangulation  or  mesh  on  with  triangles 

of  size  h,  as  shown  in  Fig.  1.  The  vertices  of  the  triangles  ^ 
are 

y(  Th _ fi 


rArjnnwi 

VArATArArAVA 

rArArArATATA 

nnnnnxi 

VArATArATATA 


VArAVArArAVA 


1  X 


Fig.  1 

called  the  nodes  of  Tj^.  By  analogy  with  the  definition  of  in 

a  ,  we  choose  as  trial  space 


Sj  »  (u:  For  each  T  « 


u 


linear  combination  of  1, 


9 


I 


m 


rx 


dt 


,  and 


ry 


ds 


b(s) 


,  u  continuous  at  the  nodes,  and 


(3.16) 


u  =  0  at  the  boundary  nodes}. 

For  the  test  space  we  choose 

=  linear  combination  of  1,  x,  and 


S-  .  =  {v:  For  each  T  6  T.  ,  v 
2  ,h  h 


y,  u  continuous  at  the  nodes,  and  u  =  0  at  the 
boundary  nodes}.  (3.17) 


In  contrast  to  the  situation  in  Example  a,  the  nodal  constraints  inposed 

on  the  functions  in  S,  .  do  not  inply  the  functions  are  continuous 

X  /  n 

and  we  have  i.e.,  the  trial  space  is  nonforming.  ^2  h  ^  ^2 


in  both  examples. 
We  now  define 


”l,h  =  “2,h  =  «l,h  =  «  L2(fi),  u 


6  Hj^(T)  for  all  T  €  Tj^}, 


lullu  “  llullu  =  ilulL  =  [[  u^dxdy  +  I  [  |Vu|^dxdy] 
”l,h  2,h  “l,h  Jfi  T€Tj^  Jt 


1/2 


and 


Bjj(u,v) 


I  f 

T€Th 


+  bu  V  )dxdy 


X  x 


y  y 


B.  is  defined  on  H,  .  ^  h,  .  and  (2.3)  holds  with  M 
h  l,h  l,h 


=  6. 


The  approximate  solution  is  then  defined  by 


"h  ^  ®l,h 


Bh{Uh»v)  =  I  fvdxdy,  Vv  €  S2 


(3.18) 


It  is  possible  to  show  that  (2.5)  holds,  i.e.,  that 


inf 

u€S 


sup  |B.  (u,v)  I  ?  y(oi,S)  >0  \/h 


(3.19) 


l,h  ^^S2,h 


where  Y(afS)  depends  on  a  and  8  but  not  otherwise  on  a(x)  and 
b(y)  .  We  have  also  shown  that  the  functions  in  ^  approximate 


the  unknown  solution  well.  In  fact 


inf  llu  -  xt! 


x«s 


l,h 


l,h 


C{a,8)h"u" 


(3.20) 


L,h 


where 


=  {u:  u  €  L2(n),  u  €  Hj^(T),  au^,  bu^  f  (T)  VT  €  T^,  u  is 
continuous  at  the  nodes  of  T^,  and  u  =  0  at  the 
boundary  nodes} 


and 


Hull 


I  (bUy)y|^}dxdY. 


It  is  easy  to  see  that  u,  au  ,  bu  €  H, (T)  implies  u  is  continuous 

X  y  J. 

on  T  and  so  the  requirement  of  continuity  at  the  nodes  makes  sense. 

Now,  combining  (3.19),  (3.20),  and  (2.9)  we  have 

||u  -  u.  ||„  <  C(a,8)h||u||  for  all  h.  (3.21) 

^  ‘^l,h  “L.h 


We  have  been  able  to  prove  the  following  regularity  result  for  (3.14) 
(or  (3.15)):  If  f  €  L2 ,  then  u  €  and 


Hull 


<  C(a,6)  HflL 

h 


(3.22) 


where 

=  {u:  u  €  H^(f}),  bu^,  €  Hj^(n)} 

and 

"“''h  “  I  I  (bUy)yl  ^}dxdy. 

Because  it  is  immediate  that  HuH„  =  ’lull„  for  any  u  €  ,  from 

”L,h  “l  ^ 

(3.21)  and  (3.22)  we  get 

Theorem  b.  The  error  u  -  u^  satisfies 

l|u  -  u,  H„  <  C(a,6)  IlfH^  Vh.  (3.23) 

"  "l,h  ^2 

(3.23)  shows  that  we  have  first  order  convergence  in  the  "energy  norm," 
with  an  estimate  that  is  uniform  with  respect  to  the  c)a«?«;  of  coeffi¬ 
cients  satisfying  0  <  a  <  a(x),b(y)  "5  . 

If  , .  .  .  form  the  usual  basis  for  S2  defined  by 

«.(z.)  =  6..,  for  all  interior  nodes  z  ,...,z  of  T.  ,  and 


form  the  basis  for  S,  .  defined  by  (f.Cz.)  =  6..,  for  a] 

X  #  n  13  X  j 

nodes  z ^ ,  then  the  stiffness  matrix  of  (3.18)  is  given  by 


where  a,  and  b.  are  the  (one  dimensional)  piecewise  hai 
— n  ~n 

averages  of  a(x)  and  b(y) ,  respectively.  The  result  in 
completely  analogous  to  that  expressed  in  (3.11)  and  (3.12) 
a  .  We  exphasize  that  although  we  are  treating  a  two  dimer 
problem,  the  harmonic  averages  are  one  dimensional. 

As  we  have  seen,  trial  spaces  with  good  approximation  e 
are  required  in  an  accurate  method.  It  is  thus  natural  to 
trial  functions.  We  have,  in  fact,  been  able  to  show  that 
defined  in  (3.16),  is  an  optimal  approximation  space  in  th€ 


sup 

“'"L.h 

'""-L.h-l 


inf 

X(S 


u 


l,h 


sup  inf  l|u  -  xHu 


tor  all 
nodes  of 
(cf.  [8]) 


.  with  dim  S 
N  L,h 

Therefore 

relative  to  the  norms 


„  =  dim  S,  V.  =  N  =  the  numbe 
N  l,h 

is  optimal  in  the  sence  of 

II  •  i|„  (the  norm  we  ai 

”l,h 


the  error  in)  and 

bounded  by  Hfll.  ) 
^2 


H 


,  ,  (a  2—  derivative  norm  in  which  i 

Xi  ^ 

We  have  thus  chosen  an  "ideal"  approx: 


subsoace  for  our  problem.  We  have  chosen  the  test  space  £ 

defined  in  (3.17),  so  that  B,  (u,v)  can  be  calculated  fror 

h 

data  for  ^  h  that  the  inf-sup  cc 

(2.5)  holds.  These  features  of  the  choices  for  the  trial  « 
space  lead  to  the  accuracy  of  the  method.  We  further  note 
choice  of  S2  led  to  an  easily  computed  stiffness  matri: 
the  precise  statement  of  the  accuracy  and  robustness  of  th« 
We  note  that  (3.25)  shows  (in  the  case  a  =  b  =  1)  that 
S-  .  of  continuous  piecewise  linear  functions  is  an  optims 

2  y  n 

matino  subspace  in  the  case  when  the  differential  operator 
Laolacian.  Furthermore,  in  an  asymptotic  sense,  S2  is  < 
oroblems  with  smooth  coefficients. 


c.  ^  Singular  Two  Point  Boundary  Value  Problem 
Consider  the  problem 


n 


& 


-(/x  u') '  =  f (X) ,  0  <  X  <  1 
u(0)  =  u(l)  =  0. 


(3.26) 


=  H  =  (u:  I  /x  (u')^dx  <  u(0)  =  u(l)  =  0] 

0 

llull^  =  llull^  =  llull^  =1  /x  (u’)^dx. 


The  variational  formulation  of  (3.26)  is  given  by 


u  (  H 

rl 

■ 

✓x 

Ev 

J 

0 

fc 

Given  a 

'  v*dx 


-r 


fvdx  Vv  €  H. 


jj  =  {u:  I  /x  (u')^dx  <  »  for  j  *  l,...,n, 


u(0)  =  u{l)  =  0} 


lullj  =  ilullj  =  llullj  + 

1 ,  h  2  2 


I 

3=1  ^ 


/x  (u')^dx. 


We  define  the  trial  and  test  spaces  by 

S,  w  =  {u:  For  each  j,  u  =  a  linear  combination  of  1  and 

l,h  ^  J 

/x,  u  continuous  at  the  nodes,  and  u(0)  =  u(l)  =  0} 

and  S-  .  =  C®  piecewise  linear  functions,  as  in  Example  a  . 

4  f  n 

The  approximate  solution  '*h  ^  ®1  h  then  characterized  by 


B  (u  ,v)  =  fvdx  Vv  €  S. 
h  h 


where 

B 


^(u,v)  =  I  u’v'dx. 

I  -I 


IV»V1ST»3 


One  can  show  that  (2.3)  and  (2.5)  hold  with  M  and  y  independent  of 
h.  Furthermore,  one  can  show  that 


inf  llu  -  xll„  ^  Cl.  a,. 

X«l,h  0 


Combining  these  facts  with  (2.9)  leads  immediately  to 

Theorem  c.  The  error  u  -  u.  satisfies 
-  h 

llu  -  u,  IL  <  Ch(f  -1^-1—  dx)^/^  Vh. 


(3.27) 


We  emphasize  that  (3.27)  holds  for  an  arbitrary  mesh.  It  is  not 
necessary  to  refine  the  mesh  in  the  neighborhood  of  the  singular  point 
0.  The  rate  of  convergence  in  (3.27)  is  the  highest  possible.  This 
follows  from  the  theory  of  N-widths.  Problems  similar  to  (3.26)  have 
recently  been  considered  in  [6,11]. 

d.  A  Boundary  Value  Problem  in  a  Domain  with  a  Corner 
Consider  the  model  problem 

-Au  =  f (x,y) ,  (x,y)  €  Q 

u  »  0,  (x,y)  €  r  =  an  (3.28) 

where  n  is  a  polygonal  domain  with  a  convex  angle  (see  Fig.  2) ,  which 


Fig.  2 

we  assume  is  placed  at  the  origin,  and  f  f  Cq(:^).  The  solution  u 
of  (3.28)  is  singular  at  (0,0)  and  as  a  consequence  the  standard 
finite  element  method  with  piecewise  linear  approximating  functions 


and  a  quasiuni form  mesh  gives  an  inaccurate  approximate  solution.  In 
fact 

-  %\(n)  " 

with  p  <  1  depending  on  the  angle  B.  (If  S  «  it,  then  p  =  1.) 

Appropriate  refinement  of  the  mesh  in  the  neighborhood  of  the  concave 

-1/2 

angle  leads  to  the  optimal  convergence  rate  0(N  ),  where  N  is 

the  number  of  degrees  of  freedom  (the  dimension  of  the  trial  and  test 
space) .  A  different  way  to  achieve  accuracy  for  problems  with  corners 
is  to  augment  the  trial  space  with  singular  functions  (which  will  not 
have  local  supports).  See,  for  example  [5,7].  We  will  outline  an 
approach  that  preserves  the  local  nature  of  the  finite  element  method 
bv  selecting  special  trial  functions.  The  resulting  method  will  have 
the  highest  possible  rate  of  convergence,  namely  0(h) . 

Toward  this  end  suppose  is  a  quasi-uniform  triangulation  of 

Q,  let 

S,  w  =  {u:  For  each  T  €  T.  ,  u  is  a  linear  combination  of  1, 

1 ,  h  n 

T 

r**  cos  a9 ,  and  r°*  sin  a0,  where  a  -  ^/B  and  (r,0) 
are  polar  coordinates  of  (x,y) ,  u  is  continuous  at 
the  nodes ,  and  u  =  0  at  the  boundary  nodes  > , 

and  let  S.,  .  be  the  usual  C**  piecewise  linear  functions  relative 
A  f  n 

to  T.  that  vanish  on  T.  The  triangulation  will  be  required  to 
n 

satisfy  certain  technical  restrictions  in  addition  to  being  quasi¬ 
uniform.  We  describe  these  now.  Let  0  <  t  and  0  <  h  <<  1  be 
given.  Choose  K  sufficiently  large,  to  be  precise,  choose  K  >  K(B,y), 
where  K(a,T)  is  an  explicitly  known  expression.  Triangulate  the 
portion  of  n  near  the  concave  angle  with  several  isoceles  triangles 
such  as  shown  in  Fig.  3.  Triangulate  the  rest  of  1  with  triangles 


Fig.  3 


of  size  h  and  minimal  angle  >  t  as  in  Fig.  4.  The  bilinear  form 
Bj^,  the  approximate  solution  u^^,  and  the  norm  '  •  are  defined 


Fig.  4 


in  the  obvious  way.  Then  we  can  verify  (2.3)  and  (2.5)  and  use  (2.9) 
to  prove 


Theorem  d.  The  error  u  -  u.  satisfies 


lu  -  u. 


’ch|log  h|^/^,  B  =  2 
Ch,  6  <  2  . 


(3.29) 


We  thus  have  a  first  order  estimate  for  the  "energy"  norm  with  a 
nearly  uniform  mesh. 


e.  An  Interface  Problem 
Consider  the  problem 

-div  (a  grad  u)  =  f(x,y),  (x,y)  €  n 
u  ■  0  on  r  =  3(2 

trhere  (2  is  a  convex  polygon  and 


(3.30) 


where  aj^  and  82  are  positive  constants  and  ?2j^  c  (2,  3(2^^  is  smooth, 
and  Q2  *  ^  "  ^1*  first  derivatives  of  u  have  jumps  across  the 

interface  3(2]^  and  thus  the  standard  finite  element  method  is  inaccu¬ 
rate  unless  the  mesh  is  chosen  so  that  the  interface  lies  on  edges  of 
triangles  in  the  trinagulation.  We  will  now  show  that  maximal  accuracy 
can  be  achieved  without  aligning  the  mesh  with  the  interface  if  we 
modify  the  trial  functions  properly. 

Let  be  a  quasi-uniform  triangulation  or  mesh  on  Q.  We 
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m 


describe  a  function 

U  €  s,  .  . 

1  ,h 

i)  For  T  €  T^, 

if  T  n 

3Q.  =  (j> ,  u  1 

is  linear 

i  't 

if  T  n 

3f2,  ^  (J),  u| 

is  as  follows: 

Choose  0  ^  n  T,  Denote  the  tangent  and  normal 

to  by  t  and  n,  respectively.  Then 

are  linear, 


a)  =  u 


TflfJ, 


and  U2  =  u 


6)  Uj^(Q)  =  U2(Q)  , 


3u.  3u- 

TC^'O'  =  ^<0>- 


3u, 


3u. 


ii)  u  is  continuous  at  the  nodes  of  T^. 


Note  that  i)<S)  is  a  requirement  that  the  trial  functions  approximately 


3u 


model  the  interface  condition  aj^ 


3u 


n. 


=  a. 


~3h  ~  “2  3n  satisfied  by  the 

exact  solution.  For  test  space  we  choose  the  usual  C”  piecewise 
linear  functions.  Having  defined  these  spaces  we  then  define  the 
aonroximate  solution  by 


%  ^  ®l,h 


B. 


(u,v)  H  T  (  aVu  •  Vvdxdy  + 

^  Jrnn- 


Tfin. 


aVu  •  Vudxdy)  = 


fj,  yy  '  S2,h' 


(3.31) 


For  the  method  defined  in  (3.31)  we  have  proved,  with  the  aid  of 
(2.9),  the  estimate 


!u  -  u.  lb  <  Chllf!'  Vh 
”  "l,h  ^2 


(3.32) 


where 


Iv'!- 

H 


l,h 


f  2 

!  \f  e 

u 


j  vudxdy  +  ( 


Vvl^dxdy  + 


Vv I ^dxdy) 


Tna, 


(3.32)  is  a  first  order  estimate  for  the  "energy"  norm  of  the  error. 
We  emphasize  that  there  is  no  relationship  between  the  mesh  and  the 


jAIiA 


interface.  The  method  proposed  here  would  provide  an  alternative  to 
methods  in  which  the  interface  is  modeled  by  a  mesh  line,  as  in  19]. 


4 .  An  Approach  to  the  Design  of  Finite  Element  Methods  for  Problems 
with  Rough  Input  Data 

Our  treatment  o^  the  Examples  in  Section  3  suggests  an  approach 
to  the  design  of  finite  element  methods  for  problems  with  rough  input 
data.  In  this  section  we  outline  the  steps  in  this  approach. 

a)  Choose  the  space  of  right  hand  sides  or  source  terms  f  to 


be  considered.  In  all  of  our  example  we  chose  f  (  L. 


Example  3.c  where  we  considered  f*s  satisfying 


Jn  /Z 


except  in 
dx  <  ». 


0  /x 

b)  Find  the  space  of  solutions  corresponding  to  the  space  of 
right-hand  sides  under  consideration.  This  will  involve  a  regularity 
result.  For  example,  in  Example  3.b  ,  the  space  H.  is  the  space  of 

U 

solutions  corresponding  to  f's  in  L^.  Often  regularity  results 
are  available  for  problems  with  smooth  data,  but  are  not  available 
in  sufficient  generality  for  problems  with  rough  data. 

c)  Choose  the  mesh  dependent  bilinear  forms  and  spaces  (norms) 


V  ''i,h<' 


'H 


l,h 


«2,h< 


'H 


2,h 


) .  This  choice  is  usually  very 


natural,  following  directly  from  the  basic  variational  formulation, 
considered  triangle  by  triangle. 

d)  Select  trial  spaces  which  have  good  (optimal)  approximation 
pronerties.  This  is  the  major  problem.  Usually  such  spaces  are 
closely  related  to  local  solutions  of  the  equation  under  consideration. 
In  many  situations  the  proper  choice  leads  to  non-conforming  functions. 
Inter-triangle  continuity  is  enforces  only  at  the  nodes.  The  approxi¬ 
mation  properties  of  the  trial  functions  is  directly  tied  to  the  space 
of  solutions  (see  b) ) .  The  problem  of  selecting  optimal  trial  functions 
is  often  not  simple,  in  practice  one  would  li)ce  to  find  a  trial  space 
that  performs  almost  as  well  as  the  optimal  trial  space  but  which  is 
easily  implemented. 

e)  Select  a  test  space  so  as  to  ensure  the  inf-sup  condition 
is  satisfied  with  constant  y  that  is  not  too  small  and  so  that  the 
stiffness  matrix  can  be  calculated.  In  contrast  to  the  trial  space, 
the  test  space  is  chosen  to  be  conforming. 

We  have  illustrated  these  steps  on  the  relatively  simple  examples 
discussed  in  Section  3.  We  restricted  our  attention  to  first  order 
methods  in  which  the  maximal  rate  of  convergence  was  0(h).  As  mention¬ 
ed  in  the  introduction,  the  ideas  in  the  paoer  can  be  generalized  in 


Tg  ^  ^ v» '.  ^  V^  '.■%  rv  ri  *"^3:^  t.^  *.  '.^  ■' 
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various  directions. 
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